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ABSTRACT 



We use a numerical code called PAKALMPI to compute synthetic spectra of 
the solar emission in quiet conditions at millimeter, sub-millimeter and infrared 
wavelengths. PAKALMPI solves the radiative transfer equation, with Non Local 
Thermodynamic Equilibrium (NLTE), in a three dimensional geometry using a 
multiprocessor environment. The code is able to use three opacity functions: 
classical bremsstrahlung, H~ and inverse bremsstrahlung. In this work we have 
computed and compared two synthetic spectra, one in the common way: using 
bremsstrahlung opacity function and considering a fully ionized atmosphere; and 
a new one considering bremsstrahlung, inverse bremsstrahlung and H~ opacity 
functions in NLTE. We analyzed in detail the local behavior of the low atmo- 
spheric emission at 17, 212, and 405 GHz (frequencies used by the Nobeyama 
Radio Heliograph and the Solar Submillimeter Telescope). We found that the 
H~ is the major emission mechanism at low altitudes (below 500 km) and that 
at higher altitudes the classical bremsstrahlung becomes the major mechanism 
of emission. However the brightness temperature remains unalterable. Finally, 
we found that the inverse bremsstrahlung process is not important for the radio 
emission at these heights. 

Subject headings: radiative transfer equation, solar radio emission, numerical model 
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Introduction 



The radio emission from astronomical sources may be a powerful tool to study the 
physical conditions of the source and the medium between it and the observer. In particular 
at centimetric, millimetric and sub-millimetric wavelengths the radio emission from the 
solar atmosphere carries out valuable information from the physical conditions of the 
chromosphere and the transition region. 

The solar emission at high radio frequencies can be explained by two major solar 
regimes: the flaring case, where non-thermal emission from accelerated electrons and 
thermal emission from hot and cold plasmas are the main processes; and the quie t Sun 



1970|), for 



regime, where the emission is attributed to thermal bremsstrahlung (iKundiJ 
which the magnetic field is not directly important. Notice that the magnetic field structures 
the matter in the solar atmosphere, therefore it will influence the plasma density structure, 
as instance, in the case of chromospheric spicules. 

In the early 1930's, the first attempts were made to explain the solar chromosphere 



flCillie fc Menzel 



19351 ). But it was only after the eclipse observations of 1952, that 



the existence of a t wo-component m odel in temperature was established (hot and cold 



temperature model. 



Giovanelli 



1949|). 



The classical treatment of the solar chromosphere in quiet regime is based on the 



assum ption that the atmosphere is static and in hydrostatic equilibrium ( jVernazza et al 



19731 ). Also, the radio and UV spectra must be in agreement, therefore in several works 



have been studied the relationship between 
millimeter and sub-raillime t er wavelengths ( 



;he UV line emission and the continuum at 



Chiuderi Drago et al. 



1983; 



Chiuderi fc Chiuderi Drago 



Kuznetsova 



Landi fc Chiuderi Drago 



2004 



de La Luz et al. 



1978 



2003; 



20081). 



Ahmad fc Kundu 



Loukitcheva et al. 



1981 



2004 
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By assuming an emission mechanis m, the observed radio-spectrum can be correlated 



with the atmospheric height of the source (jVernazza et aL 



2008 



De la Luz et al.l 



1976 



Ayres 



1989 



Avrett &: Loeser 



20101). Generally, the height associated to a specific emitting frequency 
is the height of the layer where the atmosphere becomes optically thick. However, the 
transition between optically thin (tj, ^ 1) and optically thick regimes (r^ ^ 1) is quite 
ambiguous. 

Theoretically, the solar emission can be modeled by solving the radiative transfer 
equation taking into account the temperature and density of the solar atmosphere, 
and using the required opacity functions according with the relevant photon-particle 
interactions. In particular, from millimeter to infrared wavelengths range, the free-f ree 



inter a ction (bremsstra 



1985 



Raulin &: Pacini 



ilung) is the main mechanism used to reproduce solar spectra (JDulk 



20051 ). and, also, a fully ionized atmosphere is commonly assumed. 



With these assumptions both the time and the complexity of the computations are greatly 
reduced, specially for the quiet Sun emission case. 

It is clear that more photon-particle interactions take place in the solar atmosphere, 
i n particular, the H~ co ntinuum provides the dominant radiation at photospheric level 



(^Athay fc Thomas 



19611 ) . and therefore H must be considered as an important source 



of sub-millimeter and infrared emission at low chromospheric level. We expect that the 



interactions between H and free electrons (inverse bremmstrahlung, 
1980l ) is important at this wavelenght range also . 



Golovinskii &: Zon 



Therefore, in order to reproduce a more realistic solar spectra, it is necessary to 
compute explicitly the density of the different ion species and feed with this information 
the relevant opacity functions, using Non Local Thermodynamical Equilibrium (NLTE). In 
this way, we have constructed a synthetic spectrum of the solar quiescent radio emission at 
short wavelengths. 
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As a first step, we compute the ion and electron density at any layer of tlie atmospliere 
using as inputs, atomic models and pre-calculated density and temperature atmosplieric 
profiles for the quiet Sun regime. Then, we solve the radiative transfer equation for 
the required ray paths . We use PAKALMPI, an updated version of the PAKAL code 



( JDe la Luz et al. 



20101 ) to compute a synthetic spectrum of the solar emission from 
millimeter to infrared wavelengths. PAKALMPI solves the radiative transfer equation, in 
NLTE conditions, in a three dimensional geometry using a multiprocessor environment. 

The computed synthetic spectra can be used to study the detailed behavior of the 
emission and absorption throughout the solar atmosphere. In particular, in this work we 
are interested in the amount of local emission at different atmospheric heights for the 
millimeter - infrared wavelength range. We focus at 17, 212, and 405 GH z, because these 



are t he observing frequencies of the Nobeyama Rad io Heliograph (NoRH, 



Kaufmann et al 



Nakajima et al 



20081 ). For this, we 



19951 ) and the Solar Submillimeter Telescope (SST, 
use NLTE computations for Hydrogen and H'^ ion, taking into account three opacity 
functions (bremsstrahlung, H and inverse bremsstrahlung) for twenty ion species in 
different ionization stages. 

The goal of this work is to compare the common synthetic spectra obtained by using 
the bremsstrahlung opacity function only with fully ionized atmosphere, and a more realistic 
spectra calculated by using NLTE and three opacity functions, computed with high spatial 
resolution (of around 1 km) . 

The three opacity functions, the computation of ionization stages of twenty ion species 
considering H, H~, and rie in NLTE are included in our new model called "Celestun". 
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1. Opacity Functions 



In order to include the contribution of different ion species, we use the following 
(complete) opacity function which resumes the slightly different expressions for 



the bremsstrahlung 



Rybicki &: Lightman 



opaci t y reported in the hterature (JKurucz 



1979 



Bulk 



1986 



Zheleznyakov 



19961): 



1985 : 



where: e is the electron charge; c is the velocity of light; m the electron mass; k the 
Boltzmann constant; Z^ is the charge of the ion species ^; T is the electron temperature; v 
the frequency; ng the electron density; n^ is the density of the ion species ^; and gff{Z^, T, u) 
is the average free-free Gaunt factor which takes into account the initial and final quantum 
energy states of the total momentum. 



by 



The rigorous quantum me chanical expression for the Gaunt factor has been derived 



Sommerfeld fc Maud (119351 ) but the solution is difficult to compute and therefore, 



many approximations have been developed in order to make the computation easier (see 



Brussaard fc van de Hulst 



( 1l962l ) for a review of average Gaunt factor in the bremsstrahlung 



opacity function). 



Following iRybicki fc Lightman! (119861 ) we compute the average free-free Gaunt factor 



for chromospheric conditions as: 



9ffiT,Z, 



^In 
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1, 
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X <l,y < g{x) 
x<l,y > 9ix),y < h{x) 
X < l,y > h{x) 
X > l,y < 1 
X > l,y > I 



(2) 
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where 



X 



kT 
Z'^Ry 



(3) 



and Ry = 13.6 eV is the Rydberg constant and 



fcT' 



g{x 



9999 / 99 



-X 



99 
103 



' W V102 

h{x) = 1001 - lOOOx 



(4) 

(5) 
(6) 



Under chromospheric conditions it is necessary to include the contribution of the H ion 



in the opacity expression. Th is contribution is effective through three mechanisms ( jJohn 



1988 



Golovinskii fc Zon 



1980l ): 



The Wildt's photo-detachment mechanism (jWildtlll939l ) 



hv + H' ^ H + e 



The neutral interaction (jZheleznyakovl Il996l ) 



hv + e' + H ^ H + e 



The bremsstrahlung for H or inverse bremsstrahlung fJGolovinskii fc Zorull980l ) 



hu + e~ + H~ ^ H~ + e~ 



In the milli meter - sub-millimeter wavele ngth range, the photo-detachment mechanism 



is not relevant ( lAlexander Sz Ferguson 



19941 ). whereas the neutral interaction is very 



important. The corresponding opacity function is 



^ra ~ ^u "el^H, 



where fc*°* is the total absorption coefficient pubhshed in iJohnI ( 119881 ). Pg is the electronic 
pressure and uh is the Hydrogen density. Finally, in the wavelength r ange of interest, the 



inver se bremsstrahlung opacity is the same than the positive ion case (JGolovinskii &: Zon 
1980l . we use bremsstrahlung with rii = uh-)- 



Therefore, the total opacity used to compute the synthetic spectrum is 



^ff = ^// + '^ni 



"ff ■ 



2. Density Profiles 



As input parameters PAKALMPI needs the Hydrogen density, the temperature, and 



the departure coefficients (see Sec. 12.31) . With this infor mation, PAKA 



Vernazza et al. 



MPI computes the 



(119731 ). but using a 



ionic and electronic densities following the derivation of 

higher number of energy levels for each ion, i.e the Saha equation is solved not only for the 

base state. 



We define the ionization contribution function [Z function) by: 

H F 



(7) 



where 



^ = {H,He,Li,Na,K,...,E}. 



are the atoms considered in the model; A^^ is the ionization level; n^^^ is the density of the 
atom C, in the ionization level N^; and F is the highest energy level for each ion. In order to 
calculate the Z function we need an atomic model and the relative density or abundance of 
each atom. Although the relative density may depend of the atmospheric height, in a ffist 
approximation, we assume a constant relative density for all the atmospheric layers. 
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2.1. Electron Density 



In IVernazza et al.l ( 119731 ) we found the formulas for the electron density which takes 
into account NLTE computations and uses the fundamental energy state for the ions. 
However, using our definition of ionization contribution function (Eq. [7]), we have: 

-(1 - Zd) + v/(l - Zdf + MinnZ) 



rie. 



2d 



where uh is the Hydrogen density, Z is the ionization contribution function, T the 
temperature and d the NLTE contribution function. 

2.2. Ion Density 



To compute n^^^ , we solve the matrix for the classical Saha equation: 



log 






-0.1761 -log(Pe)+lo 



UN. * 



5040 



where 



UNf 



J29i- 



AT exp 



X^,Ni:,k 

kT 



(9) 



(10) 



PAKALMPl can use a single model of atoms, but is able to manage different atomic 
densities in each layer. Even more, it is possible to change the met alicity in each layer o: 



the atmosphere. In this work, we use the atomic model published by 



Vernazza et al. 



(119731). 



To calculate the density of H we use the following expression (IVernazza et al. 



19731): 



Uh- = 1.0345 X 10-^^6j:^-nenH/T-=^/2g8762/T_ 
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2.3. NLTE Contribution Function 



The NLTE contribution function is: 



d = bMT) (i + E ^) ' 



yN 



where the term ^^=2 — << 1 fIVernazza et al. 



19731), so that, 



d = hij{T), 



(11) 



(12) 



where 



^(T) 



h^ 



3/2 



ewiXn/kT), 



27nnkT ^ 
and the depar ture coefficient f or the Hydrogen in th e fundamental state (6i) defined in 



(13) 



Menzell (119371 ) and modified by 



Vernazza et al. 



(119731 ) is: 



h 



rii/nl 



where the ratio nl/nl is given by the Saha equation in thermodynamical equilibrium. 

As 6i is difficult to compute for each specific atmospheric layer and due to the fact 
that there are only few published atmospheric models reporting this parameter, we have 
developed an algorithm to obtain an approximation of the departure coefficient based on 
any published bi coefficients. 

We found that different theoretical models of the chromosphere are very close to each 
other when plotted in the dynamic space (density versus temperature, the analysis in the 
dynamic space is important because in this space, the problem of the ambiguity between 
height and temperature is avoided, this is, when the temperature proffie has two or more 
associated heights for a single temperature point). Therefore, it should be easy to get new 
approximated fei parameters from published departure coefficients. If we have two models: 
a published model (model A) of the atmosphere with: temperature T^(/i), density if^(/i). 



- 11 - 

and departure coefficient parameter bi{h), where h is the height over the photosphere; and 
if we have a new model (model B) with temperature Tb(/i) and density Hsih); we can 
approximate bi from model A as follows: 

1. At a given height over the photosphere {hs) we take the temperature and density 
from the model B 

Pb = {TB{hB),HB{hB)). 

2. With Pb we found the closer point jTa = {TA{hA),HA{hA)) to the model A (in the 
dynamic space). This point has an associated height (/i^) over the photosphere in the 
model A. 

3. We take the points that enclose Ha 

[/lo, hA, hi] 
in the same model A. 

4. With ho and hi we can evaluate hi in the model A in order to interpolate hA'- 

blihB) ^ biihA) = MhltiM^^hA - ho) + biiK). (14) 

hi -ho 

By construction, pa and pb have similar physical conditions. Therefore, the computed 
departure parameter 6i(/?-b) can be used in the model B with high degree of confidence. 
This process is applied to all points in model B to obtain the required approximation for 
the departure parameters. 
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3. Radiative Transfer Equation 

The local emission, or the emission in a given atmospheric layer, can be represented in 
the following way 

^ local ^abs ' ^emi \^^) 



where 



is the local absorption and 



eabs = Ioexp{-Ti, 



ocal ) 



(16) 



Slocali^ - exp{-Tiocal)) 



(17) 



is the local emissivity. 

In this case /q is the incoming emission from the background (lower adjacent 
atmospheric layer), Tiocai is the local optical depth and Siocai is the local source function. 

We use an update d version of PAKAL code to solve the radiative transfer equation 



( JDe la Luz et al. 



2010l ) along a ray path using an iterative method to compute Iq and /, 



local- 



The Cabs and eemi are useful to analyze the behavior of the local emissivity of the 
atmospheric layers. While eabs is related to the level of opacity of a given layer, eemi 
indicates the capacity of this layer to produce radiation. 

With this definition of the radiative transfer equation, we can study in detail the local 
emission and absorption processes. 

4. Algorithm 



The code PAKALMPI performs the following steps to compute the emission I^: 
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• Reads the temperature, Hydrogen density, and the atomic models for a given layer 
from the atmospheric pre-calculated model (VALC, C07, etc), 

• Solves the matrix of the Saha equation (Eq. [9]). The solution of the matrix is a set of 
density of ions for each atomic species considered in this layer. 

• Finds the "Z" contribution of electrons in LTE (Eq. [Tj). 

• Computes the approximation to the 61 parameter using the temperature and density 
of the specific layer. 

• Using 61 and Z, computes the initial values of rie (Eq. [8]) and n^. 

• Computes the hh- and recalculates the ng and uhi (this step is necessary since the 
computation of uh- needs the electron and neutral Hydrogen densities). 

These steps are repeated recursively until rie converges. The iterative process is necessary 
because we are searching for the equilibrium of the physical system. 

Once He has converged, the results for all ions are saved and the code gets the 
temperature. Hydrogen density, and abundances for the next layer. In this way, the code 
produces tables with the height profiles of the electron and ion densities. 

When all the atmospheric layers have been calculated, PAKALMPI computes the 
ray trajectory (from the Sun to the observer) and the intersection of this ray path with 
the radial models (ions, temperature and electronic density) from the computed tables. 
Finally, PAKALMPI solves the radiative transfer equation using all the considered opacity 
functions. 

The process is repeated (for each ray trajectory) to create 2D images at several 
frequencies. In this way, PAKALMPI improves the construction of 2D images considering 
detailed 3D emitting structures. 
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5. Chromospheric Spectra 



As inputs to our code , we use two complem entary atmospheric models, the 



chromospheric model C07 (jAvrett &: Loeser 



20081 ) which includes temperature and 



hydrogen density but does not include the departure coefficients, and the VALC model 



Vernazza et al. 



198ll ) in order to get these departure coefficients. 



The (atmospheric) plasma behavior can be easily visualized by plotting the plasma 
density versus the temperature, this is the dynamic space. Therefore, in Figure [T] we have 
plotted the neutral Hydrogen density as a function of the temperature for VALC and C07 
models. We note that the differences between both models are minimal. This gave us 
confidence to use both: the COT atmospheric model in order to get Hydrogen density and 
temperature; and the departure coefficients for the fundamental state of H {hi) and H~ 
{hu-) which are interpolated from VALC model using our approximation method (Sec. 
12. 3p . The result is seen in Figure Ej where we have plotted the original departure coefficients 
(VALC) with a dashed line and our interpolation (for COT) with a continuous line. The 
main differences correspond to the region of the temperature minimum and the transition 
region, and are caused by temperature differences between both models. 

Our code output gives us detailed information of the radiation parameters. As instance, 
in Figure |3]we show the computed optical depth at IT (top), 212 (middle) and 405 (bottom) 
GHz as a function of radial height. We made two simulations, the ffist one using only 
bremsstrahlung opacity (dashed black lines); and the second one using our model Celestun, 
in this case, the yellow lines represent the inverse bremsstrahlung opacity; the H~ opacity 
is shown by red lines; the blue lines correspond to the classical bremsstrahlung opacity; 
and the black continuous lines represent the full opacity function from Celestun. The 
introduction of H~ and inverse bremsstrahlung opacities does not change the height where 
the atmosphere becomes optically thick (r ~ 1). However, it is clear that below 500 km. 
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the H process is the major source of opacity in the solar atmosphere. 

The detailed local behavior of the emission at 17, 212 and 405 GHz is seen in Figure IH 
where we have plotted the brightness temperature (black lines) and the input atmospheric 
radial temperature (yellow lines) as a function of the atmospheric height. The continuous 
and dashed lines represent, respectively, the cases where all, and only bremsstrahlung 
opacities were considered. The contribution, to the local brightness, of the absorption eahs 
(Eq. [16]) is plotted in blue lines, whereas the emission egmi (Eq- [IZ]) is plotted in red lines. 

We note that there are two separated regions where clearly the emission has a significant 
contribution: the first one is close to the photosphere, and it is interesting to note that the 
height of this region increases when the H~ opacity function is considered; The importance 
of the second enhanced emission region decreases with the frequency, and the H~ opacity 
does not produce any differences in this region. The second enhanced region and the H~ 
opacity do not affect the final brightness temperature, this is why these structures were 
unnoticed until now. The contribution of the inverse bremsstrahlung is negligible. 

6. Synthetic Spectra 

It is interesting to compare the output of our code using different input atmospheric 
models. As instance, in Figure |5] we show observations per formed at different phases 



of the solar cycle and collected by 



Loukitcheva et al 



(120041 ): diamonds corresponds to 



close-to-minimum, the square symbols to intermediate, and triangles to close-to-maximum 
phases of the solar cycle. In the same figure, we have plotted the synthetic solar spectra 
computed using all (continuous lines) and only bremsstrahlun g (dashed lines) opa cities and 



the following input models: VALC in yellow lines; CAIUS05 (Selhorst et al. 



20051 ) in blue 



lines; C7 in black lines; and again CAIUS05 but using their electronic density and physical 



-16- 



assumptions (bremsstrahlung opacity and fully ionized atmosphere) in red lines. 

There are interesting differences in the input atmospheric models, e.g., the altitude 
and the length of the chromospheric temperature minimum, as well as the altitude and the 
morphology of the chromospheric upper limit (i. e. the beginning of the transition zone 
and corona). The VALC model has the temperature minimum closer to the photosphere 
and is the warmest. This model has a thin chromosphere and defines a "plateau zone", 
where the chromosphere temperature reach 20,000 K. This plate au is responsible of the 



high brightness temperature at low frequencies (see Figure 151 and ide La Luz et al.ll2008l ). 



The CAIUS05 model has similar altitude for the temperature minimum but is hotter and is 
shorter compared to the VALC model. The chromospheric temperature grows exponentially 
up to coronal values. The CAIUS05 model has the largest chromosphere. The C7 model 
has the lowest temperature minimum, and then a constant temperature until the rise due 
to the transition region and the corona. This model eliminates the plateau zone using 
ambipolar diffusion. 

In order to facilitate the analysis, we consider separately three regions of the spectrum: 
from 2 to 60 GHz (low frequency); from 60 to 400 GHz (intermediate frequency); and from 
400 GHz to 3 THz (high frequency). 

At low frequency range, the temperature differences of the atmospheric models 
produce significant differences in the spectra (regardless of the opacity function), this 
emitting region is located at the upper chromosphere ~ 1500 km over the photosphere. For 
intermediate frequencies, the simulations agree rather well with the observations, however, 
the observations have large dispersion. Finally, at high frequency range we note that there 
is an apparent discontinuity at 400 GHz. In this region the model that better fits the 
observations is CAIUS05. On the other hand, both C07 and VALC models predict higher 
brightness temperatures than the observed. 
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These apparent minimal differences between input atmospheric models, produce 
significant differences in the resultant synthetic spectra computed from our radiative model. 

7. Discussion and Conclusions 

Using different atmospheric models and opacity functions, we have computed the 
optical depth, r, at low atmospheric altitudes. Our results show that if we use either, only 
bremsstrahlung or Celestun (i.e. if, H~ and inverse Bremsstrahlung) opacities, the height 
over the photosphere where the overlying atmosphere becomes optically thick (r ?^ 1) 
remains unalterable. However, if we take into account the local emissivity (Eq. [T71) we show 
that the height where the main emission is generated, changes meaningfully. 

For example, at 17 GHz the height over the photosphere where the atmosphere starts 
to become optically thin is around 1500 km, but the height where r = 1 is around 2050 km. 
There are local emissivity structures below 1500 km, although they remain invisible to the 
observer because the atmosphere is optically thick at these heights. When we considered 
bremsstrahlung opacity only, the local emissivity at 17 GHz changes in a region between 500 
and 850 km over the photosphere (the atmosphere becomes optically thin). This structure 
is not present when the Celestun model is considered. 

At 212 and 405 GHz the behavior of the local emissivity changes noticeably, since 
at these frequencies two regions of emissivity are present: A lower one, closer to the 
photosphere shows deep differences when comparing between both models. The height of 
the region where the emission is generated is higher by 200 km when the Celestun model is 
used. The uppermost region is located between 800 and 1500 km above the photosphere. In 
this region there is no difference between the two models. As seen in Figure HI the second 
layer of emission is below the r = 1 height and is similar for both models. Even more, after 
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1000 km, the brightness temperature remains unalterable. This is due to the fact that the 
temperature of the atmospheric model is the upper limit for the brightness temperature, 
and remains constant between 1000 and 2000 km above the photosphere, the brightness 
temperature remains unalterable. These facts mask the local emission processes which are 
unnoticed by the observer. 

We found that the dominant emission mechanism changes with altitude: H~ is 
dominant below 500 km whereas bremsstrahlung dominates at upper heights, although the 
final brightness temperature remains unchanged. 

Our results show that the inclusion of the H~ and inverse bremsstrahlung mechanisms 
does not modify the final brightness temperature. But the H opacity function has a 
significant impact in the local emissivity and absorbency and it must be considered when 
studying micro structures in the low chromosphere. Therefore Celestun code is a valuable 
tool to perform this kind of studies. The inverse bremsstrahlung has no impact on the 
brightness temperature, neither on the local emissivity nor on the absorbency processes. 

We note that there are more emission mechanisms act ing at the low solar atmosphere, 



as instance, molecular emission has been reported for C O (ISolanki et al. 



199- 



Faurobert fc Arnaud 



and also in 



20021). The 



the 515.60 to 516.20 nm spectral range (C2 and MgH, 
emitting layer has been situated below the 1100 km, making evident that these molecules 
are present in a thin layer inside the temperature minimum region, between the photosphere 
and the low chromosphere. However, observations in visible wavelengths, have shown that 
the molecular emission at the limb brightening region is negligible, corresponding to 0.2% 
of the continuum intensity at the center of the solar disk. Nevertheless, a detailed study of 
the contribution of the molecular rotational spectrum at these wavelengths is necessary. 

Although our simulations are related to quiescent solar situations, the next natural 
steps will be to investigate and compare with flaring cases. Therefore, such researches 
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could in principle be used to better understand the thermal structure of the flaring solar 



atmosphere, that is the response of the atmosphere to 



existing semi-empirical flare models like 



f are en ergy release. From the 



Machado et al. 



fll980f ) and 



Mauas et al 



fll990[ ) we 



ew 



should grossly expect decreasing radio fluxes as a function of increasing wavelengths in the 
range from centimeter to far-infrared wavelengths. 

We will pay special attention to the density structure of the flaring loops (due to 
its importance in the computed opacity functions). It is clear that new simulations, 
observations and their comparison are needed to better investigate the structure and 
dynamics of the low layers of the solar atmosphere during flares. 
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Fig. 1. — Temperature vs Hydrogen density for VALC model (solid line) and C07 Model 
(dashed line) in the dynamic space. 
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Fig. 2. — Departure coefficients for VALC (daslied line) and C07 (continue line) approxi- 
mated model taking into account the bi parameter published in the VALC model. We note 
that for both models 6i parameter follow the radial temperature. The main differences are 
seen above 1500 km over the photosphere. 
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Fig. 3. — Optical Depth at 17, 212, and 405 GHz using the C07 input atmospheric model. In 
blue we show the contribution of the bremsstrahlung opacity, red H~ and yellow the Inverse 
bremsstrahlung. The continuous line shows the sum of all the contributions and the dashed 
line represents only the bremsstrahlung opacity function. 
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Fig. 4. — Local absorption [tabs-, blue) and local emissivity (eemi, red) for 17, 212, and 405 
GHz (top to bottom) using the C07 input atmospheric model. In yellow we show the radial 
temperature and in black the brightness temperature. 
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Fig. 5. — Sy nthetic solar spect r a and observations at different stages of solar activity cycle, 



collected by 



Loukitcheva et al. 



(J2004J ). The square points are intermediate, diamonds close- 



to-minimum, and triangles close-to-maximum of solar activity cycle. The solar spectra are 
shown for bremsstrahlung opacity (dashed line) and Celestun model (continuous line) for 
four input atmospheric models: C7 (black), SEL05 (blue), VALC (yellow) and SEL05 but 
using their electronic density (red). 



